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1 . Formulation  of  the  Problem 

Consider  a communication  system  consisting  of  r 
stations,  each  of  which  must  be  able  to  communicate  with  all 
the  others.  The  communication  is  conducted  via  an  intermediate 
satellite.  Since  each  station  has,  realistically,  a finite 
capacity  to  handle  messages  simultaneously  in  progress,  and 
since  the  satellite  itself  has  limited  capacity,  the  system 
will  sometimes  be  congested,  and  a message  applying  for  trans- 
mission will  be  blocked,  i.e.  effectively  given  a busy  signal. 
We  wish  to  calculate  the  probability  that  a message  will  be 
blocked,  or  delayed.  The  reader  familiar  with  telephone 
system  congestion  theory,  see  Syski  [I960],  or  Cooper  [1972], 
will  recognize  this  as  a more  complicated  version  of  the 
situation  for  which  the  "Erlang  B"  formula--a  truncated 
Poisson — holds.  Essentially  we  assume  that  blocked  calls  are 
lost.  For  a model  describing  the  more  realistic  "re-try" 
situation,  see  Gaver  and  Lehoczky  [1976]. 


D.  P.  Gaver  and  J.  P.  Lehoczky  acknowledge  the  research 
support  of  the  Office  of  Naval  Research  at  the  Naval  Post- 
graduate School.  We  are  grateful  for  helpful  discussions  with 
Dr.  Martin  Fischer  of  the  Defense  Communications  Agency. 


Specifically,  we  assume  that  station  i (1  <_  i <_  r) 

has  c.  channels,  and  the  satellite  has  c„  channels.  We 
1 s 

assume  that  a message  initiated  at  station  i and  intended  for 
station  j (i  to  j for  short)  requires  a free  channel  at  i, 
one  at  the  satellite,  and  one  at  j before  transmission  may 
begin.  If  no  channel  is  available  at  one  of  these  three 
locations,  blockage  occurs.  We  assume  the  satellite  offers 
direct  access,  thus  if  any  channel  is  available  in  the  satellite 
a user  will  not  be  blocked  at  that  point.  Furthermore  we 
assume  that  a channel  in  use  is  not  available  to  any  other 
user.  That  is,  there  is  no  possibility  of  simultaneous 
transmission  by  another  user  on  an  occupied  channel  and  con- 
sequent message  spoilage.  The  possibility  of  message  destruc- 
tion apparently  exists  for  some  existing  satellite  communication 
systems;  see  Kleinrock  [1975],  and  Gaver  and  Lehoczky  [1977]. 
Finally,  direct  access  structure  is  apparently  not  yet  avail- 
able in  practice,  according  to  our  information.  Our  study 
pertains  to  conceptual  systems. 
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A 


2 . Analysis 

Suppose  that  attempts  to  transmit  messages  from  i 

to  j arrive  according  to  a homogeneous  Poisson  process  with 

rate  A^,  the  rate  of  message  termination,  when  the  call  is 

from  i to  j,  is  p^,  and  holding  times  are  independently 

exponential . Let  p . . = A . ./u ■ • , and  n • • = p..  + p . • . Let 
p Hij  13  13  13  il  31 

X^j(t)  be  the  number  of  messages  or  calls  in  progress  from 
i to  3 at  time  t.  It  is  clear  from  our  formulation  that 
= {X^ j (t) , 1 £ i < j £ n)  is  a multivariate  Markov  process 
in  continuous  time. 


Steady-State  Solution 

Note  that  satisfies  inequality  constraints  C 
which  occur  because  the  channel  capacity  is  limited  at  the 
various  stations. 


C:X12(t)  + X13(t)  + •••  + X L ( t ) + X21(t)  +•••+  Xrl(t)  < Clf 


and  in  general, 


I X (t)  + l X (t)  < c. 

j*i  13  )*L  31  1 

l l X (t)  < c 
i=l  ^ S 


1 < i < n 


(2.1) 


3 


It  the  constraints  were  not  present  (c^  = ■ , cg  - ) then 

X ^ ^ ( t ) is  an  infinite  server,  Poi sson  arrival  queueing  process 
(termed  M/M  ■)  for  every  station  pair  i,  j,  and  the  stationary 
distribution  is  then  Poisson i 


lim  PiX  . ( t ) = n..}  = e 


n . ■ 0,1,2,.. 


(2.2) 


Fur t hermore , the  number  of  calls  in  proqress  between  all  pairs 
of  stations  are  independent.  It  may  even  be  stated  that  the 
above,  (2.2),  is  true  for  the  arbitrary  service  time  situation. 
If  the  constraints  are  seldom  binding,  that  is  if  blockinq  is 
a rare  event,  then  (2.2)  provides  a useful  approximation. 

In  the  case  that  the  constraints  are  imposed,  the 
above  result  is  also  very  nearly  true,  as  is  seen  from  the 
fol low l nq : 


Result.  The  stationary  joint  distribution  of  is  Poisson 

constrained  to  the  region  C.  That  is. 


lim  P(X  =n)  lim  P{  X,  _( t ) =n X , (t)=n  ,} 

t t t 12  12'  r , r- 1 r,r-l 


n . . n . 

. e.1]  / -n.  . n.11 

. w = n n e 'J  -LL/V  i,n  e ^ -L2_ 
n n . . ! / * . , . n . ! 

i^l  n ' ■'  i/l  13 


(2.3) 


if  n r f,  and  is  zero  otherwise.  That  is,  if 
13 


f 


nl2+  nl  3+  ’ * * + "ir  + n21  +”*+  nr  1 £ C1 


etc. , as  in  (2.1) 


Discussion.  In  order  to  justify  the  solution  (2.3)  we  consider 
the  balance  equations  for  the  steady  state  probabilities  it. 
These  have  the  following  form  at  state  values  away  from  the 
boundaries,  the  latter  being  defined  by  the  constraint  set  C. 


tt  (n^ » • • • ' * * * ,n2r'  * " * ' *^r  1 ' " " ” * ^r , r~  1 ^ 


* |X12  +‘”+  Xr.r-1  * 


= l h (. . . ,n  +1, . . . ) ( (n . .+l)u . . ) 

i/j  13  3 3 


+ l n ( . . . . n . . - 1 

i?j  13 


(2.4) 


I 


These  equations  state  that  the  rate  of  departure  from  the 
state  n equals  the  rate  at  which  that  state  is  entered. 
Actually,  there  is  local  balance ; if  n(n^_.)  denotes  marginal 
distribution  of  calls  in  progress  between  i and  j,  then  in 
the  unconstrained  case  we  can  see  that  local  balance  holds: 
if 

it  (n)  = II  n (n . . ) (2.5) 

i*j  3 
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Then,  termwise  in  (2.4)  for  all  i / j, 


n ( ... ,n  \ \ . + n . . ii  . . 

i]  i]  lj 


= ,i(...,n..+l,...)(n  . + l)M  + ir(...,n..-l,...)A..  (2.7) 

i]  ij  il  il 


since 


‘pij  (> 


i] 


n . . ! 
il 


A.  . = e"'  ^ 
il 


n . + 1 

<P,J  1J 


(ni  .+1) ! (ni j 


+1) v . . 
il 


il 


n . . - 1 
, il 


il  _LL 


n . . ! 
il 


n . . vi  . 
ilil 


il  _2J_ 


diij-lM 


A . 


il 


This  shows  that  the  product  from  solution  (2.5)  holds  for 
n strictly  within  C.  Now  suppose  n is  a boundary  point. 


This  means  that  some  transition  rates  which  were  A. . > 0 

il 


(2.8a) 


(2.8b) 


in  the  unconstrained  case  must  be  equal  to  zero,  in  order 
to  keep  the  X process  within  f,  i.e.  on  the  left  hand  side 
of  the  balance  equations  (2.4)  these  terms  now  involve  zeros 


fttaVV'' 


**  V 


"V 


lllgjfr;  u,|  . ..  . v-  — - 


No.  of  Copies 


0 


Now  inside  C the  balance  equations  (2.10)  are  satisfied 
by  the  product  form 


(n12'n21) 


12 

'P12  P12 
n12  ‘ 


^ < 0 1 2 ) T,^n2i^ 


21 

21  f 12 
n2 1 ' 


(2.11) 


Now  suppose  n^9  + n9^  - min (c^ ,c2  >c^) , i.e.  is  on  the  boundary 
of  i.  Then  • ^ 7 + ■ , ^ must  be  set  equal  to  zero.  But, 
correspondingly  v (n12  + l ,n., } ) - n(n  ,n21+l)  = 0.  By  local 
balance,  the  product  form  solution  continues  to  hold  on  the 
boundary.  Write  for  n12  + n2  - boundary  point 


(0  + 0 


f n 1 2 u 1 2 + n21p21 


(2.12) 


and  cancel  off  common  factors;  the  balance  is  obvious.  It 
is  only  necessary  to  normalize  the  product  form  over  the 
constraint  region,  as  dictated  by  (2.3). 
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No.  of  Copies 


3 . Blocking  Probabilities 

The  probability  that  a call  originating  at  station  i 

is  blocked,  essentially  receiving  a busy  signal,  is  calculated 

in  principle  from  (2.3).  It  is  convenient  to  define 

Y^U)  = X^lt)  + X ^ ± ( t ) for  1 < i,  j£r*  i / j and 

Y^(t)  = 0.  Here  Y^^  (t)  represents  the  total  number  of 

calls  in  progress  between  stations  i and  j.  In  steady  state 

Y„  are  independent  Poisson  random  variables  with  parameter 

n..  = p . . + P-.,  constrained  by  C : 
ij  13  31 ' 1 


J. 

y Y . . < c . , l<i<r 
j=l  ^ - 1 - “ 


^ l l Y.  . < 
l L-  L.  1]  - 


(3.1) 


i 3 


3 - °s 


Now  observe  that  a call  from  i to  j can  be  blocked 
in  three  ways: 

1)  At  the  originating  station,  if  Station  i is  full. 

This  is  event  E.  = {J^  . Y..  = c.}. 

2)  If  the  satellite  channels  are  full,  the  event 

E = {T.V.  y. . = 2c  ). 
s ‘-l^j  ij  s 

3)  If  the  destination  station.  Station  j,  is  full.  This 
is  event  E ^ . 

The  probability  an  i to  j or  j to  i transmission  is  blocked 
somewhere  is  given  by 
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P(i  to  j or  j to  i transmission  is  blocked) 

= P(E.  U E . U E ) 

13s 

= P(E.  ) + P ( E . ) + P(E  ) - P(E.  n E . ) 

1 j s 13 

- p (e . n e ) - P(E  n E ) + p (e . n e^  n e ) . 0.2) 

1 o Jo  1 J t) 

Each  of  the  above  probabilities  can  be  represented  in  terms 
of  the  Y. . random  variables.  The  value  of  each  of  these 
probabilities  can  be  easily  found  by  summing  terms  of  the 
form  (2.3),  the  steady  state  distribution,  over  a boundary 
portion  of  C.  For  example 


while  other  terms  in  (3.2)  can  be  computed  by  changing  the 
numerator  to  reflect  a change  in  the  boundary  conditions. 

It  is  clear  that  the  calculation  of  each  of  the 
terms  in  (3.2)  is  in  principle  straightforward  as  it  involves 
merely  the  calculation  of  a well-defined  ratio.  Unfortunately, 
the  problem  may  be  nearly  computationally  infeasible  if  the 
c.'s,  cg,  and  k are  large.  For  example  if  c^  = cg, 

1 < i < r,  then  C includes 


cs  + 9 

(k) 

1 2 


distinct  points.  If  cg  = 50,  then  for  k = 3,4,  and  5 this 
quantity  is  2.3426  x 10^,  3.2468  x 10^,  and  7.5394  x 10^ 
respectively.  Many  interesting  cases  are  essentially  com- 
putationally infeasible. 

Computer  programs  have  been  written  for  the  cases  of 

r = 3 and  4 (r  = 2 can  be  done  with  the  Erlang  B formula) . 

It  is  possible  to  reduce  the  computations  necessary  in  (3.3) 

as  follows.  Let  c . = min(c.,...  c ,c  ) . The  denominator 

min  i.  it  s 

(and  numerator)  can  be  rewritten  as 


rkt.  c 


l l n n = I s 

n=0  ytC  k < Z yk«.  n=l  n 


Z E ykfn 
k < Z K 


c ■ 
min 


= y s + y S . 

1 _ n l _ n 


n=0 


n=c  . +1 
min 


; 3 . 4) 


Now  using  the  multinomial  theorem 


cmin 

l s 

Ln  n 
n=0 


Cmin 

l 

n=0 


*kS. 

Ou  n 

l n n 

yeC  k < l ykZ' 

^ ^ yk£=n 
k < Z 


C . 

min 

i 

n=0 


l 


n n 


yklt 

’kfc 


EE  y.  0=n  k < Z *kr 
k < £ 


c . ( E E nk») 

^in  k < Z k* 


n=0 


n ! 


(3.5) 
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The  last  term  is  simply  computed.  This  observation  removes 


c . + 

min 


(k) 

2 ' 


(k> 
' 2 ; 


can  reduce  the  computations  required  substantially.  Neverthe- 
less, for  interesting  values  of  k,  cg,  and  c^,  1 < i < r, 
the  number  of  terms  needed  to  be  computed  may  render  the 
method  to  be  infeasible.  Research  directed  toward  finding  a 
tractable  approximation  useful  for  large  networks  is  presently 
under  way. 
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Numerical  Results  for  Three  and  Four  Stations 


We  now  present  a few  numerical  results  that  have  been 
obtained  for  the  situation  in  which  three  or  four  ground 
stations  communicate  via  satellite.  The  computer  programs 
used  for  obtaining  these  numbers  is  available  upon  request. 

It  calculates  the  probabilities  using  enumeration  of  the  multi- 
nominal  terms.  Three  stations  require  a relatively  small  number 
of  computations.  For  the  case  of  four  stations,  the  reduction 
(3.5)  is  utilized. 

We  are  interested  in  cases  where  the  blocking  prob- 
abilities are  small,  say  less  than  .10.  We  wish  to  see  if  in 

such  circumstances  probability  of  blocking  (P(E.  1)  E.  II  E )) 

1 J s 

can  be  estimated  assuming  independence.  Specifically,  we  wish 

to  determine  if  P(E.  U E.  u E ) can  be  approximated  by 

1 3 s 

1 - P(E^)  P(Ej)  P(Eg) . If  such  an  approximation  is  reasonable, 
it  reduces  the  amount  of  computation  required  in  the  problem. 

In  looking  over  the  following  tables,  it  appears  that  this 
approximation  is  usefully  accurate,  especially  for  the  cases 
of  small  (less  than  .1)  block  probability. 


Case  3, 


r = 4,  c1  = c2  = c3  = 10,  nij  = 2.0 


fs 

Exact 

Approxima 

10 

.3019 

.3023 

12 

.1991 

.2038 

14 

.1226 

.1335 

16 

. 0805 

.0916 

18 

.0683 

.0733 

20 

.0674 

.0692 

Case  4. 


r = 4,  c.  = 


C2  c3 


10,  r^j  = 3.0 


c 

s 

Exact 

Approxim 

10 

.493 

.494 

12 

.402 

.408 

14 

.319 

. 339 

16 

.254 

.286 

18 

.221 

.242 

20 

.216 

.223 

15 
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